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In the present work we use the large-A^c approximation to investigate quark matter described 
by the SU(2) Nambu-Jona-Lasinio model subject to a strong magnetic field. The Landau levels 
are filled in such a way that usual kinks appear in the effective mass and other related quantities. 
/3-equilibrium is also considered and the macroscopic properties of a magnetar described by this 
quark matter is obtained. Our study shows that the magnetar masses and radii are larger if the 
magnetic field increases but only very large fields (> 10^^ G) affect the EoS in a non negligible way. 
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I. INTRODUCTION 

In 1979, telescopes in spacecrafts and astronomers 
around the world detected the emission of very intense 
gamma and X rays. The sources of these rays were first 
called soft gamma repeaters (SGR) and later identified 
as possible remnants of supernova explosions, the tragic 
death of very massive stars. If this remnant is a neutron 
star that spins very rapidly, an intense magnetic field 
is formed; if it spins slower the magnetic field is strong 
but not as much as in the first case. The ordinary neu- 
tron stars, also known as pulsars, bear a magnetic field 
of the order of 10^^ — 10^^ G. The neutron stars with 
very strong magnetic fields of the order of 10^'* — 10^^ 
G are known as magnetars and they are believed to be 
the sources of the intense gamma and X rays detected in 
1979. Most of the time the magnetar remains inactive, 
but the strong magnetic field causes the solid crust to 
break into small pieces. The crustquake leaves a fireball 
which cools down and emits X rays from its surface until 
it evaporates completely [l|, Q ■ 

In Ref. Q the equation of state used to describe neu- 
tron stars in a strong magnetic field is obtained from 
a field theoretical approach. Two relativistic models 
are used: one normally called non-linear Walecka model 
(NLWM) and the other one with a derivative coupling 
between mesons and baryons. The importance of includ- 
ing anomalous magnetic moments (AMM) is discussed. 
A more recent work Q analyses the importance of the 
scalar-isovector 5 mesons in the EoS that describes mag- 
netars. In [5] density dependent hadronic models are 
used with the same purpose. In all three works the AMM 
of the electrons were not considered because they were 
shown to cancel out if properly introduced [H. In two 
papers [1,11 the AMM of the muons were also taken into 
account. All above mentioned papers refer to neutron 
stars composed of hadrons and leptons. These stars are 
called hadronic stars. 

In the stellar modeling, the structure of the star de- 
pends on the assumed equation of state built with appro- 



priate models. The true ground state of matter remains 
a source of speculation. In conventional models, hadrons 
are assumed to be the true ground state of the strong 
interaction. However, it has been argued [1, 0, S, M, [l3| 
that strange quark matter (SQM) is the true ground state 
of all matter. This hypothesis is known as the Bodmer- 
Witten conjecture. Hence, the interior of neutron stars 
should be composed predominantly of u, d, s quarks, plus 
leptons to ensure charge neutrality. Pulsars described 
by matter composed of SQM are often called strange 
stars. However, as the strangeness content depends on 
the model used to describe the quark matter, we prefer 
to describe any model in which the interior involves de- 
confined quarks (not bound in hyperons) as quark stars 
[Ti'|. Apart from the differences in the EoS, an important 
distinction between quark stars and conventional neutron 
stars is that the quark stars are self-bound by the strong 
interaction, whereas neutron stars are bound by gravity. 
This allows a quark star to rotate faster than would be 
possible for a neutron star. Moreover, some authors have 
argued that quark stars should be bare [13, Ell , in the 
sense that any crust would either not form or would be 
destroyed during the supernova explosion. The charac- 
teristics of the radiation from hot, bare strange stars have 
been identified [ij] and the electron- pos itron pairs that 
can be emitted from bare quark stars [l5j require the ex- 
istence of a surface layer of electrons tied to the star by 
a strong electric field. Two of the most common models 
used to describe quark matter are the MIT bag model 16 1 
and the Nambu-Jona-Lasinio model (NJL) [131 • In 11 j 



it was shown that while the electron chemical potential of 
a quark star described by the MIT bag model is very low 
(less than 20 MeV), the NJL model gives much higher 
values reaching 100 MeV inside the star, accounting for 
the necessary electric field that explains the emission of 
electron-positron pairs. 



An important point to 
stability of quark matter 
Two different possibilities 
found in relation to quark 
stars: the unpaired phase 
in the literature on stranj 
locked phase (CFL) [F 
the Fermi surface to form 



be investigated refers to the 
in the interior of quark stars, 
for the MIT bag model can be 
matter in the interior of quark 
Toj. which is widely favored 
ge stars, and the color-flavor- 
which allows the quarks near 
Cooper pairs which condense 
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and break the color gauge symmetry [l9| . At sufficiently 
high densities the favored phase is the CFL. In the 
stability of quark matter described by the two-flavor NJL 
model and subject to an external magnetic field was in- 
vestigated. It was shown that the stability depends on 
the strength of this field and on model parameters. 
In Ref. 2l| the EoS for magnetized quark stars was 



where the quark sector is described by the Nambu-Jona- 
Lasinio model 



described with the help of the MIT bag model. AMM for 
the quarks were properly taken into account. 

The NJL model in a magnetic field was considered in 
[l^l and it was shown that the magnetic field sponta- 
neously breaks chiral symmetry. The formalism used in 
the above mentioned study was based on a previous cal- 
culation performed in Ref. [23j . 

The scope of the present work is to study ud quark 
matter in a magnetic field within the NJL model, with 
or without the requirement of /3-equilibrium. We focus 
our work on the SU(2) version of the model. Moreover, 
we remark that the present work may also be relevant 
regarding the physics of non-central heavy ion collisions 
such as the ones performed at RHIC and LHC-CERN 
which can provide a possible signature for the presence of 
CP-odd domains in the presumably formed quark-gluon 
plasma phase 124|. In this particular case one reaches 
magnetic fields of about 10^^ G or B ~ Qni^/e (m^r rep- 
resenting the pion mass and e the fundamental electric 
charge) . 

This work is organized as follows: In section II we 
set up a lagrangian density adequate to describe two fla- 
vor quark matter, in (3 equilibrium, in the presence of 
an external magnetic field. This allows the derivations 
to be carried out in a uniform fashion from a common 
lagrangian density making the evaluations more trans- 
parent and pedagogically more interesting to the reader. 
Using functional techniques we derive the effective poten- 
tial, which represents Landau's free energy density, using 
the standard large-iVc approximation. In section III the 
relevant expressions for quark matter EoS subject to a 
magnetic field are displayed. Specific cases, for the sym- 
metric version, are discussed in section IV A for finite 
density and no magnetic field, IV B for zero density in 
a magnetic field and IV C for the general case of finite 
density matter in a magnetic field. In section V, details 
of matter in /3- equilibrium are given and in section VI 
our results are shown and discussed. In section VII the 
more important conclusions are drawn. 



II. GENERAL FORMALISM 

In order to consider (two flavor) quark stars in (3 equi- 
librium with strong magnetic fields one may define the 
following lagrangian density 



(2) 



where a summation over the quark flavors, f = u,d is 
implied while qf represents the quark electric charge. 

Note that we have used rric = ~ md as representing 
the current masses. 

The leptonic sector is given by 



(3) 



1 /? FA^" 



(1) 



where Z = e, /i. One recognizes this sector as being repre- 
sented by the usual QED type of lagrangian density. As 
usual, and F^^^, — d^Ai, — di,A^^ are used to account 
for the external magnetic field. Then, since we are in- 
terested in a static and constant magnetic field in the z 
direction, A^ = 5^2XiB. 

Regarding the actual evaluations, £/ and £/ bear two 
fundamental differences. Due to the quadratic fermionic 
interaction, the former is non renormalizable in 3-1-1 di- 
mensions (G has dimensions of eV~^ ), meaning that 
eventual divergences cannot be eliminated by a con- 
sistent redefinition of the original model parameters 
(fields, masses, and couplings). The renormalizability 
issue arises during the evaluation of momentum integrals 
which represent Feynman loops and, in the process, one 
usually employs regularization prescriptions (e.g. dimen- 
sional regularization, sharp cut-off, etc) which are formal 
ways to isolate divergences. However, the procedure in- 
troduces arbitrary parameters with dimensions of energy 
which do not appear in the original lagrangian density. 
Later, when (unlike the 3-1-1 NJL model) the theory is 
renormalizable, one may choose any value for the arbi- 
trary energy scale and the original parameters run with 
it as dictated by the renormalization group. Within the 
NJL model a sharp cut off (A) is preferred and since the 
model is nonrenormalizable one gives up the very high 
energy scales fixing A to a value related to the physical 
spectrum under investigation. This strategy turns the 
3-1-1 NJL model into an effective model while A is treated 
as a parameter. The experimental values of quantities 
such as the pion mass (m^r) and the pion decay constant 
(/tt) are used to fix both, G and A. 

A second important issue regards the fact that, when 
TOc , the quark propagator brings unwanted infra 
red divergences meaning that the evaluations have to be 
carried out in a nonperturhative fashion. Moreover, very 
often physical quantities (such as the self energy) appear 
as powers of the dimensionless quantity Gh? which is 
greater than the unity spoiling any possibility of success 
via standard perturbative evaluations. 

As far as analytic nonperturhative evaluations are con- 
cerned, one can consider one loop contributions dressed 
up by a fermionic propagator whose effective mass (M) is 
determined in a self consistent way. This approximation 
is known under different names, e.g., Hartree, large- iVc 
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or mean field approximations (MFA). The leptonic sec- 
tor, on the other hand, is described by a QED type of 
lagrangian density which is renormahzable and which, 
in principle, could be treated in a perturbative fashion. 
Here, we treat the complete Cpf in a consistent way by 
evaluating the thcrmodynamical potential related to Cpf 
up to one loop. It is interesting to remark that, in prac- 
tice, one does not have to consider the full Lpf since 
both sectors have similar polynomial structures, apart 
from the four quark interaction term which, as discussed, 
makes the theory nonrenormalizable and renders pertur- 
bative calculations useless. So, concerning the evaluation 
of the equation of state, our strategy is the following. We 
use quantum field methods, in the imaginary time formal- 
ism, to evaluate the effective potential, or Landau's free 
energy density (.F/), for the quark sector in the large- 
Nc approximation. After summing over the Matsubara's 
frequencies we will regularize the divergent (three) mo- 
mentum integrals by using a sharp non covariant cut off, 
A. When evaluated at its minimum, the effective po- 
tential is equivalent to the thcrmodynamical potential, 
VLf = —Pf=£f — TS — fj,fpf where P/ represents the 
pressure, £f the energy density, T the temperature, S 
the entropy density, and /i/ the chemical potential (a 
sum over repeated indices is implied). The quark den- 
sity, p f , and the baryonic density, p b , are simply related 
by p/ = 3pB- For the present study, just the zero temper- 
ature case is important and, as a consequence, the term 
with the entropy vanishes. Equivalent results for the lep- 
tonic sector, within the same approximation, can be triv- 
ially obtained by performing the replacements G — > 0, 
rric mi, qf ^ qi and Nc ^ \- Finally, this procedure 
allows us to obtain the EoS from the full pressure for (3 
stable dense quark matter in the presence of a magnetic 
ficld(B) 

Ppfipf.PuB) = Pf{pf,B)\M + Pi{puB)U,+^ , (4) 

where our notation means that Pf is evaluated in terms 
of the quark effective mass, Af, which is determined in a 
(nonperturbative) self consistent way while Pi is evalu- 
ated at the leptonic bare mass, to/. The term 5^/2 arises 
due to the electromagnetic term Fi^^^F^^ /A in the original 
lagrangian density. In order to normalize our results we 
shall require that the pressure vanishes at zero chemical 
potentials by defining 

PpfM^^f,PuB)=Ppf{pf,puB)-Ppf{0,Q,B) . (5) 

Note that the above normalization prescription, which 
washes away the i?^/2 term, is not unique and one could 
as well require that the pressure vanishes at zero chemical 
potential and zero magnetic field in which case the i?^/2 
term survives. Although we use the former prescription 
for most of the time we will also consider the latter in 
section VI, when stellar matter is discussed. Throughout 
this paper we consider the following set of parameters 
dl: A = 587.9 MeV , to^ = 5.6MeV, to^ = 0.511 MeV , 
= 105.66 MeV, and GA^ = 2.44. 



III. EOS FOR QUARK MATTER AT FINITE 
DENSITY IN A MAGNETIC FIELD 

To obtain the effective potential (or Landau free en- 
ergy density) for the quarks, Tf, it is convenient to con- 
sider the bosonized version of the NJL which is easily 
obtained by introducing auxiliary fields (cr, 7?) through 
a Hubbard-Stratonovich type of transformation. Here, 
J-f is evaluated using the large- iVc approximation which 
is equivalent to the mean field approximation (MFA). 
Then, to introduce the auxiliary bosonic fields and to 
render the theory more suitable to apply the large-A'c 
approximation it is convenient to use G X/{2Nc) for- 
mally treating Nc as a large number which is set to its 
relevant value, Nc = 3, at the end of the evaluations. To 
retrieve some well known results related to chiral symme- 
try breaking (CSB), in the absence of a magnetic field, 
let us set B = for the moment ^. One then has 

^V^/l* 1^) V'/-V'/('7' + «75T7r)'0/--^(o-^ + 7f^) , (6) 

where a' = a+rric. Note that the introduction of the aux- 
iliary (or background) fields does not change the physics 
since they do not propagate and their Euler-Lagrange 
equations of motion trivially lead to a = —(X/Nc)'ipfipf = 
—2Gipfipf and tt = — 2GzV'/75tV'/. However, their intro- 
duction simplifies the selection of the relevant contribu- 
tions at any order in 1/Nc [1^. Basically, each closed 
quark loop contributes with a factor of Nc while each in- 
ternal bosonic line brings a factor of 1/A^c- The effective 
potential (or Landau's free energy density), Tf, is defined 
as the classical potential which appears (at the tree level) 
in the original plus radiative (quantum) corrections. As 
it is well known, J-f is particularly useful in the study of 
symmetry breaking/restoration since a symmetry which 
is observed at the classical level may be broken by quan- 
tum corrections with the appearance of a non vanishing 
order parameter. The effective potential, which is also 
the generating functional all IPI contributions with zero 
external momenta [2^ . can be readily obtained by in- 
tegrating over the fermionic fields. Within the large- A^c 
approximation this procedure results in |i27i] 

= + « J (^tr In y - (a' + tj^fTr)] . 

(7) 

The first term on the right hand side is the classical 
(tree) contribution while the second accounts for a radia- 
tive (loop) contribution of order- A'c. Let us define some 
important physical quantities by quickly reviewing how 
CSB arises within the NJL model. For this let us take 
the chiral limit (toc = 0) . Equation ([7|) can be written as 



As we shall see in the sequel one may easily incorporate con- 
tributions for T, fj, and B by modifying some of the relevant 
Feynman rules 
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In [-p2 + 



(8) 



where we have defined x = -v/c^ + tt^. Since the potential 
is symmetric we can set tt = considering only the a 
direction. Minimizing with respect to a leads to the well 
known gap equation 



da 



= 



(9) 



which, as expected, gives the self consistent relation 



1 



(2^) 



(10) 



Chiral symmetry is broken when the true ground state 
lies in (a) ^ 0. Therefore (a) is the order parameter 
which signals CSB which comes as no surprise since, by 
the equations of motion, (a) = —2G{ipf4'f) and the ex- 
istence of a non vanishing quark condensate breaks CS. 
Equation pO)) also shows that (a) is identical to the quark 
self energy within our approximation which allows us to 
write the effective quark mass as M ~ rric + {cr). Fi- 
nally, in the present work, the relevant quantity is the 
thermodymic potential, 17/, which is defined as Landau's 
free energy at its minimum — !Ff{{a)). The equation 
of state can be obtained by using Vlf = —Pf- Then, Eq. 
([5]) allows us to write the quark pressure, away from the 
chiral limit, as 



Pf = - 



[M - mcf i 



4G 



2'' 



ln[V + A/'] . (11) 



In order to obtain results valid at finite T and in the 
presence of an external magnetic field B one can use the 
following replacements, which come from dispersion re- 
lations for quarks |28i 



of 2qB due to the field in the z direction while the energy 
associated to linear motion along z is taken as a contin- 
uous. All these levels for which the values of H- lie 
between 2qBn and 2qB{n + l) now coalesce together into 
a single level characterized by n whose number is given 
by 



S 



{2^y 



dpxdpy 



SqB 
~2^ ' 



(12) 



where S is the area in the x—y plane and q stands for \qf\. 
Then, using those replacements, taking the trace, and 
summing over Matsubara's frequencies (see Appendix A) 
one obtains 



Pf 



AG 



s.nj ^ ' 



„(B)+M/]/Tl 



+ Tln[l + e-[^-(^)-^/l/^]} , (13) 



where Ep{B) = y/p^ + {2n + 1 - s)\qf\B + AP. Next, 
by analyzing the degeneracy of the lowest Landau level 
(LLL) one can define Ep^klB) = y/pl + 2k\qf\B + AP 
replacing n by fc in the sum which appears in Eq. p3p 
which, now, also runs over the degeneracy label, = 2 — 
dko- Being mainly concerned with the case T = 0, /i 7^ 0, 
and B 7^ we can take the limit T -> in Eq. (fT3)) 
arriving at (see Appendix A) 



(M- 



4G 



Tjmed 



2ti 



d 00 |.oo 7 

Y^Y.^,(\q,\B) j ^^^Ep,,{B), 



(14) 



where the contribution from the medium is 



^ pi + {2n + l- s) , with s = ±1 , n 0, 1, 



d'p T\qf\B 



i2ny 



271 



E 



dpz 



^ J (27r) 



In the above relations, = (2;/ -f- 1)ttT, with i/ = 
0, ±1, ±2, . . . representing the Matsubara frequencies for 
fermions while n represents the Landau levels (LL) and 
s represents the spin states which, at B 0, must be 
treated separately. As explained in Ref. [1^ one may 
understand the origin of those replacements, by recalling 
that, quantum mechanically, the energy associated with 
the circular motion in the x—y plane is quantized in units 



praed 



27r 



Y^Y^^MB) 



f=u k=0 



dpz 
(27r) 



[fif-Ep^B)] 



(15) 

Note that, although not explicitly written, a 9 function 
(more specifically, d{fJ-f — Ep ^)) must be considered as 
multiplying all ^/ dependent terms, such as PJ^<^'^^ ap- 
pearing in our work (see Appendix A). Simple power 
counting reveals that the last term in Eq. (fT4|) is (ultra vi- 
olet) divergent while the second term, which contains the 
in medium contributions, is finite since it has a natural 
cut off given by the Fermi momentum, p^^ p = plj — M^. 
In Appendix B we show how Eq. (|14p can acquire a 
physically more appealing form by separating the (di- 
vergent) vacuum contribution from the (finite) magnetic 
field contribution. As a byproduct, those manipulations 
also produce more elegant relations in which the infinite 
sum over Landau levels appearing in the last term of Eq. 
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(IT4)l are shuffled into Riemann-Hurwitz C functions. We 
finally get 



(Af - me)^ 
4G 



pvac pmag _^ pmec 



where the vacuum contribution reads 



(16) 



(17) 



with Ep = ■\/p2 + K'P and Nj = 2. We are now in po- 
sition to present the explicit expressions for P™'', p™"^^ 
and pp<^d.. Let us start with the vacuum which, upon 
using a sharp non covariant cut off, A, can be written as 



r}vac 



87r2 





[(A + ea)] 




M 



eAA[A' 



(18) 

where we have defined ea = VA^ + M^. The evaluations 
performed in Appendix B also give the following finite 
magnetic contribution 



pmag 



f=u 

2^ f 



1 x'-f 
- -\Xf - Xf]lnxf + ^ 



(19) 



M = TOc + 
M2 



MGNcNf 



■In 



^fAv/A^TlP 
(A + VA2+M2)2' 



M2 



2 

•■^ — ' TT 

f=u 

iln(2^)+x/-^[2x;-l]ln[x/] 



f=u k=Q 



M\qf\BNcG 



X < In 



(22) 



where we have used some (^{z, a) properties given in Ap- 
pendix B. Note that our Eq. although obtained in a 
different fashion, exactly agrees 2 with the result obtained 
in Refs. 0,1231 where the authors have used Schwinger's 
proper time formalism [s^ . In order to make this paper 
self contained let us review, in the next section, some 
standard results obtained for symmetric quark matter. 



IV. EOS FOR SYMMETRIC MATTER 



where Xf = A'P/{2\qf\B) while C{-l,Xf) = 
dC,{z, Xf)/dz\z=-i where C(z, Xf) is the Riemann-Hurwitz 
zeta function fsi'l. Finally, after integration, the medium 
contribution can be written as 



E E 

S/(fc,P)2ln 



{fif^fi}-Sf{k,By 



QflBN, 

Sf{.k,B) 



,(20) 



In this section we use Eq. (|T6)) to reproduce some of 
the most important results concerning symmetric quark 
matter at and/or B ^ Q. The result for this 

particular case can be readily obtained by setting = 
= It is worth emphasizing that, when the magnetic 
field is turned on, the u and d quark densities are not the 
same, as it is discussed next. In this sense, both quark 
chemical potentials are identical, but the densities are 
not and hence, matter is not strictly symmetric. Besides 
this, we have opted to keep the nomenclature. 



where Sf{k,B) — y/ + 2\qf\Bk. The upper Landau 
level (or the nearest integer) is defined by 



A/2 



P},F 



2\qf\B 2\qf\B- 



(21) 



A. EoS for symmetric matter at /i 7^ and 5 = 

This case has been extensively discussed in the litera- 
ture and we refer the reader to Refs. [2^, [3^ for more 
details. 

By setting B = one obtains the following relation for 
the pressure 



Finally, the term M entering the quark pressure is just 
the effective self consistent mass at finite density and in 
the presence of an external magnetic field: 



When one sets /tu = fij = fi. 
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B. EoS for symmetric matter at /i = and B 0. 



NcNf 1 




[(A + ca)] 








NcNf 1 




'{ii+PfY 


8^2 1 




[ M,, \ 



A[A2 



5/^ Q O 

- M PF 



where pj? — fJ,^ — M^. The effective mass, Af^, satisfies 



M„ = m, + ^ <i AeA - M2 In 



IJ,PF + In 



A + £A 



IJ-+PF 



(23) 



The normalized pressure is simply given by Pf csifJ-, 0) 
^/(a*,0)|m, -P/(0,0)|mo where 



P/(0,0) 



[Mq - mj2 iV.jV^ 
4G 87r 
AeA,o(A2 + ei^o)} 



2 i^^oln 



A + EA.O 



A/n 



(24) 



with eA,o = vMf+A^. The effective quark mass, AIq, 
appearing in Eq. ([M]) is just evaluated at fi — 0. 
Namely 



Mo - mc + — j^-^ < AeA.o - Mq In 



A + eA,o' 



Mo 



(25) 

Finally, at B = and /i 7^ 0, the equation of state for 
the quarks is given by 5/(/i, 0) = — P/^cff (Mj 0) + /^p(Mj 0) 
where /3(/i, 0) = dPf(fi, 0)/(i/i is the mean field result 



37r2 ^'^ '^^ 3^2 



We can also define the bag constant 



(26) 



(27) 



As emphasized in Ref. [25| one should remark that, in the 
same way as in the bag model, B describes the pressure 
difference between the trivial and non trivial vacuum, 
but its is not an input of the model, being a dynamical 
consequence of the interactions which leads to vacuum 
masses M(0, 0) 7^ TOc- Using our chosen values for G and 
A one obtains B =(181 MeV)^. Regarding CSB, one ob- 
tains the quark effective mass, M ~ 400 MeV for fi — 
and observes a first order transition for chiral symmetric 
matter at /Zc — 360MeV. A quantity of particular inter- 
est for the present work is the energy per baryon as a 
function of the density. The relations obtained above for 
the B — case will allow us to discuss, at the end of this 
section, the influence of the magnetic field regarding the 
stability of quark matter. 



In this subsection we set /i = and concentrate in the 
behavior of the quark condensate under the influence of 
an external magnetic field. One of the most remarkable 
effects of this external field regards its role so as to en- 
hance chiral symmetry breaking. This issue is related 
to the phenomenon known as magnetic catalysis which 
has been well exploited by Klimenko and collaborators 
among others [l^, . Going back to our general relation 
for the pressure, Eq. (fTB)) . and setting /i„ = = ^ = 
one obtains that the relevant pressure is given by 



P/(0,S) = - 



[M 



B ~ rUc 



AG 



Tjvac , pmag 



Mb 



(28) 



with the effective quarks mass at /j, = and B 0, AIb, 
being determined by 



Mb = Tllr 



— ^In 



(A + VA2 + M2)2 



EIVj 
- 

f=u 



«''l"l™=<=fl„r|x,|-il„(2,) 



+.,-i|2.,-ll,„.,} (») 

where now Xf = Af|/(2|g/|B) (see Eq.dH])). 
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FIG. 1: Quark effective mass as a function of B showing how 
the latter enhances CSB (magnetic catalysis). 



Fig. [T] shows the behavior of the effective quark mass 
as a function of B indicating that the former increases 
(although slightly) with the latter which stabilizes the 
condensate so that the gap equation always has a non 
trivial solution for finite B. As pointed out in Ref. [s^ 
the B field facilitates the binding by antialigning the he- 
licities of the quark and the antiquark, which are then 
bound by the NJL interaction while an electric field, E, 
has an opposite effect opposing condensate formation by 
polarizing the ipfi^f pairs. 
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C. EoS for symmetric matter at /i 7^ and B 7^ 0. 

In this subsection, still considering only the symmetric 
case, we shall review how B stabilizes quark matter due 
to magnetic catalysis, according to Ref. [20|. For this 
we can consider the general relation, Eq. ifTH)) with /Xu = 
fid — fj,, from which the quark density is easily extracted 
as 

p{f,,B) = J2 ^'^-^^^kpj(k,sf) , (30) 

f=u k=0 

where kF,f{k,Sf) = fi^ ~ Sf{k, B)"^ while the effective 
pressure is defined by 

PJ^,s{^^,B) = Pf{fjL,B)\M ~ Pf{Q,B)\Ms ■ (31) 
Then, the energy density follows as 

£pf{fi, B) = -P/3/,cff(/i, B) + mp(m, B). (32) 

In all the above equations M is given by Eq. ((22l) with 
the obvious substitution fiu = fJ-d = I^i Fig- HI 
we present some results which compare the energy per 
baryon {E/ps) as a function of the density for £? = 
and B = 2 X 10^^ G. 




1050 

1000 ^ - 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Pb/Po 

FIG. 2; Energy per baryon as a function of the baryonic 
density for symmetric quark matter. Continuous line B — 
2 X 10^" G, dashed line B = 0. 

It is shown that for i? = the curve has an absolute 
minimum at ps — 2po and that at pB =0 there is a 
local minimum corresponding to the real empty physical 
QCD-vacuum. This fact indicates that, in this case, the 
vacuum is metastable but the situation changes drasti- 
cally in the presence of an external magnetic field and 
the graph has a maximum at = so that at arbitrar- 
ily small densities stable quark droplets appear in the 
system. One then concludes that the external magnetic 
field supports the creation of stable quark matter. Later, 
when discussing our numerical results, in section VI, we 
again address some issues related to the symmetric case. 
As one can see, in figure [21 the minima obtained for both 
cases have different values. The case of symmetric mag- 
netized quark matter has a minimum at a density lower 
than the one obtained for symmetric quark matter. The 
concavity of the curve in that case is also greatest. 



Now, after reviewing some of the most important issues 
regarding how symmetric quark matter is affected by the 
presence of an external magnetic field we can incorporate 
/^-equilibrium in our results to consider stellar matter. 

V. ASYMMETRIC QUARK MATTER WITH 
/3-EQUILIBRIUM IN A MAGNETIC FIELD 

In a star with quark matter we must impose both, 
P equilibrium and charge neutrality [13]. Through out 
this paper we only consider the latest stage in the star 
evolution, when entropy is maximum and neutrinos have 
already diffused out. The neutrino chemical potential is 
then set to zero. For /3-equilibrium matter we must add 
the contribution of leptons (electrons and muons) in a 
magnetic field to the energy density and pressure. The 
relations between the chemical potentials of the different 
particles are given by 

Pd^ Pu+ Pe, Pe = Pfj^- (33) 

For charge neutrality we must impose 

Pe+ Pt, = '^('2pu- Pd)- (34) 

Now, Eqs. ([55)1 and ([M]) will have to be satisfied together 
with the self consistent relations for the quark effective 
mass during the numerical evaluations. The EoS for the 
leptonic sector is also needed. Let us now obtain this 
EoS by recalling that, as emphasized in the introduction, 
the total leptonic pressure, Pi{pi, B), is quickly recovered 
from the quark pressure , Pf{pf,B), upon performing 
obvious replacements such as / ^ and Nc = 1. Also, 
because the leptonic sector does not have the analog of 
the quartic interaction, G 0, so that when translating 
the results one takes M rric mi which lead to 

Pi ipi , B) = [pr + pr' + pr"] ™, • (35) 

Being evaluated at the same (one loop) approximation 
level all terms have exactly the same mathematical struc- 
ture as those in Pf{pf,B) (including the divergences in 
the vacuum contribution). A major difference is that all 
terms in Eq. ()14p are evaluated with the bare mi re- 
flecting the fact that an undressed lepton propagator has 
been used. At this one loop level of approximation the 
leptonic contribution is that of a free gas of relativistic 
fermions and the divergences contained in the vacuum 
contribution can be properly absorbed with a zero point 
subtraction. However, note that according to our nor- 
malization procedure we require Pi{pi, B) = at pi — 
as in the quark case which leads to the following effective 
pressure for the leptonic sector 

PiMl^i,B) - [Piipi,B) - Pi{0, B)]^^ - , (36) 

since Pi{0, B) = p^ac _^ pmag ^j^j thcsc quantities are 
written in terms of the /x;-independent mass, mi. The 
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result shows that, at the one loop level, only the following 
(finite) medium contribution has to be considered 



l=e k=0 ^ 



Slik, BY In 



si{k,B) 



(37) 



Then, the leptonic density is easily evaluated yielding 



pi{fii,B) = ^ ak-^kF,i{k,si) , (38) 

l=e k=0 

where kp^i{k, si) = \/ fJ^f ~ si{k, B)"^. Finally, the leptonic 
energy density reads 



-Pi,caOM:B) + ^iipi{^ii,B) , 



(39) 



where, again, a sum over the repeated (I) indices is im- 
plied. Finally, the total effective pressure corresponding 
to the theory described hy Cpf in the presence of a con- 
stant external magnetic field, B, is 

Pf3fMf^f,Pi,B) = Pf^^s{pj,B) + Pi^^g{fii,B) . (40) 

We now have all the ingredients to evaluate the EoS in 
/? equilibrium since 

^f3fifJ'f,f^i,B) = -Ppf,cs{fJ'f,f^i,B)+nipi+HfPf , (41) 

where the quark and leptonic densities are given by Eqs. 
(|30l) and jMl) respectively. 



VI. RESULTS AND DISCUSSION 

We start by discussing, in more detail, the effects of the 
magnetic field on the EoS of symmetric matter described 
by the usual SU(2) version of the NJL model without 
/3-equilibrium whose inclusion will be addressed after- 
wards. Whenever mentioned in the figures, Bq = 10^^ 
G. By symmetric matter we mean ud matter for which 
the chemical potential of both particles are equal. For a 
zero magnetic field this is equivalent to having symmet- 
ric matter. For finite magnetic fields due to the different 
electrical charge of both quarks and the appearance of 
Landau levels, ud is symmetric only for restricted densi- 
ties. 



In the next four figures, i.e.. Fig. [SKlSb,!!^ andHh, the 
results for B — 0.2 x 10^^ G are almost coincident with 
the results for magnetic free matter. They are kept so 
that minor differences can be seen. 

In figures [5^ and we compare the effective quark 
mass and the baryonic density as a function of the chem- 
ical potential for matter subject to different values of 
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FIG. 3: a) Effective quark mass and b) baryonic density as 
a function of the chemical potential for matter subject to 
different values of the magnetic field, po was chosen as 0.17 
fm"^ and Bo = 10^^ G. 



the magnetic field. One can see that a magnetic field 
of the order of 0.2 x 10^^ G barely affects the effective 
mass as compared with the results for ordinary matter 
(not subject to the magnetic field). Due to the Landau 
quantization, the increase of the strength of the magnetic 
field provokes a decrease of the number of the filled LL 
and the amplitude of the oscillations is more clear in the 
graphics. For each value of the magnetic field, the kink 
appearing at the smallest chemical potential corresponds 
to the case when only the first LL has been occupied. 
For B — 5 X 10^^ G matter is totally polarized for chem- 
ical potentials below 490 MeV. For the small values of 
the magnetic fields the number of filled LL is quite large 
and the effects of the quantization are less visible. For 
the larger magnetic fields the chiral symmetry restora- 
tion occurs for smaller values of the chemical potentials 
which, however, correspond to larger densities as can be 
seen from Fig. [SJ)). 

In Fig. 2^ one can see that the inclusion of the mag- 
netic field makes matter more and more bound. The 
energy per baryon E/A of magnetized quark matter de- 
scribed by the SU(2) version of the NJL model is less 
bound than nuclear matter made of iron nuclei, ^ jse^g ^ 
930 MeV, which means that quark matter is not the pref- 
erential ground state matter even in presence of the mag- 
netic fields under consideration (ssl . l36l | . 

In Fig. IIJd we plot the EoS for different values of the 
magnetic field. Once again one can see that the magnetic 
fields modify the EoS but the modifications are more sig- 
nificant when the magnetic fields reach values higher than 
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FIG. 4: a) Binding energy and b) EoS for different values of 
the magnetic field for ud matter with equal chemical poten- 
tials. Bo = 10^^ G. 



0.2 X lO^^G; as can be seen these graphics show also oscil- 
lations every time that a different number of LL is filled. 
As the magnetic field increases the number of LL for a 
given interval of energy is reduced, increasing the gap 
between them and making more visible the oscillations, 
the so-called Haas van Alphen effect fSTj. 

We have identified the densities at which pressure is 
zero, for different values of the magnetic field. This cor- 
responds to the saturation density of symmetric quark 
matter when S = 0, or to /.t„ = /i^ otherwise. In par- 
ticular, for S = 0.2 X lO^^ G, the density is 1.44 fm'^ 
and for S = 2 X 10^^ G, the density is 1.19 fm~'^, and 
for i? = 3 X 10^^ G, the density is around 2.7 fm~'^, i.e. 
the saturation density decreases as the magnetic field in- 
creases for magnetic fields lower than B = 2 x 10^^ G. 
For larger fields the opposite may occur. The minimum of 
the binding energy is the result of two contributions with 
different behavior: on one hand, due to Landau quantiza- 
tion the kinetic energy decreases with the increase of the 
magnetic field and on the other, the effective bag param- 
eter defined by the interaction terms minus the vacuum 
energy increases because the restoration of chiral sym- 
metry occurs at larger densities. The saturation density 
depends on how fast each contribution changes with den- 
sity. 

We now study quark matter in /3-equilibrium, the orig- 
inal motivation of our studies aimed at understanding the 
constitution of magnetars. In what follows, NJL SU(3) 
refers to the EoS of quark matter with u, d, s quarks 



and parameters used in [Tl| and NJL SU(2) refers to 
ud quark matter. In both cases, matter is not sub- 
ject to a magnetic field. The designation (NJL SU(2))b 
refers to ud quark matter subject to a magnetic field. 
The reader should keep in mind that here the compari- 
son with the SU(3) case has only a qualitative character 
since both, the SU(2) and SU(3), versions employ differ- 
ent parametrization sets flawing any conclusion related 
to quantitative aspects. Just as in symmetric matter, 
also for matter in /3-equilibrium, the effects of Landau 
quantization are clearly seen in the EoS (see Fig. [6^ )• 

The strange quark matter within the MIT bag has been 
discussed in [2^, [3^. Once the s-quarks are removed, 
so that only u, d matter is considered one can see that 
although the overall behavior is similar, there are clear 
differences between both models. The zero pressure den- 
sity decreases with the increase of the magnetic field for 
the MIT model while there is no clear trend for the NJL. 
These densities are given in TableHl where only u and 
d quarks are taken into account. For quark stars this 
density defines the density at the surface. We should 
remember, however that it is only a field above 10^^ G 
which has noticeable effects and, at the surface we expect 
much smaller fields. 

As for hadronic stars, the density is zero for a null 
pressure at the surface if no magnetic field is included. 
The magnetic field gives rise to an increase of the density 
in the surface In this sense, the NJL model predicts 

a more similar behavior, i.e., an increase of the density 
at zero pressure as a result of the magnetic field, whilst 
the MIT model predicts the opposite. 



B{G) 





O.lBo 


Bo 


2Bo 


3Bo 


PNJL (fm^^) 
Pm/t (fm^^) 


0.40 
0.54 


0.47 
0.47 


0.48 
0.45 


0.44 
0.36 


0.54 
0.31 



TABLE L Baryon densities of /3-equilibrium matter for which 
the pressure becomes negative, for the NJL model and the 
MIT model with the bag pressure (180 MeV)^. In both mod- 
els only u and d quarks were considered. 

The saturation density is defined by the way chiral 
symmetry is restored: a chiral symmetry restoration at 
smaller energies implies a smaller saturation density. As 
discussed before, the saturation in the NJL model is a 
balance between the attractive scalar field which satu- 
rates at chiral symmetry restoration and the kinetic term 
which gains importance once chiral symmetry restoration 
has occurred. 

The conditions of charge neutrality and /3-equilibrium 
give rise to different chemical potentials, and therefore 
different quark fractions, as seen in Fig. [5] While for B — 
/3-equilibrium matter is formed essentially by two thirds 
of d quarks and one third of u quarks and only a very 
small fraction of electrons, the presence of the magnetic 
field changes the particle fractions. Matter becomes more 
symmetric for a finite magnetic field larger than 10^^ 
G, and a larger fraction of electrons and muons occurs. 
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In Fig. [5] the lepton population is not shown so that the 
effects of the magnetic field on the quark population can 
be better noticed. 



0.6 ■•" 
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^0 



u quarks 



P/Pr 



FIG. 5: Quark population {pi/p,i = d,u in terms of density 
for B = O.lBo Bo, 2Bo and 3Bo, with Bo = 10^*' G respec- 
tively for d quarks from top to bottom and u quarks from 
bottom to top. 

In Fig. [S^ we show the EoS for /3-equilibrium quark 
stellar matter. We compare the SU{2) and SU{3) NJL 
EoS for magnetic free matter with the EoS for i? = 2 x 
10^^ G and B = 3 X lO^^ G. In the last two EoS the 
effect of the LL quantization is clearly seen. Both EoS 
are harder than the corresponding B = EoS. 

Notice that, contrary to hadronic matter in /?- 
equilibrium, where the magnetic field makes the EoS 
softer [1, the NJL model predicts a slightly harder 
EoS, as seen in Fig. [6^. 

To obtain the properties of the stars described by these 
EoS, we have added the contribution of the magnetic 
field to the pressure and to the energy density, i3^/2. As 
discussed in section II the formal consideration of this 
term is related to the prescription adopted to normalize 
the pressure. At the surface the magnetic field should 
be not larger than ~ 10^^ G and, therefore, we have 
introduced a density dependent magnetic field (soj 

B{p) = B, + B,[1~ expi-aip/poT]] , 

where po is the saturation density, Bs = 10^'"' G is the 
magnetic field at the surface, Bi is the magnetic field 
at the interior for large densities, and the parameters 
a = 5 X 10~^ and 7 = 3 were chosen in such a way that 
the field increases fast with density to its central value 
but still describes correctly the surface, namely with a 
zero pressure. 

The properties of compact stars were obtained from the 
integration of the Tolman-Oppenheimer-Volkoff equa- 
tions, using the EoS obtained with the density dependent 
magnetic field and which includes the magnetic field con- 
tribution. For i? = 3 X 10^^ G the magnetic field con- 
tribution is as large as the contribution of stellar matter. 
This is seen from the properties of the star with maxi- 
mum mass: the gravitational mass becomes larger than 
the baryonic mass, (see Table II). 



The SU(3) EoS is also included. It becomes softer af- 
ter the onset of the s quark, which occurs at a quite 
high density in this model [111- The SU(3) version of the 
NJL model is obviously more complete and the chiral 
symmetry restoration of the s quark produces a smooth 
change of declination in the EoS around 6 fm^** as can 
be seen in Fig. [SI As this energy density is very high, 
the strangeness content of the star is small. The compar- 
isons with the SU(3) NJL version has to be considered 
with care because this version of the model was fitted to 
a different set of variables, which, as we can see from Fig. 
[51 describe stellar matter (below 5 fm^^) in a different 
way from the SU(2) version. 

The results of the integration of the TOV equations are 
shown in Table II and in figure [BJs. We do not take into 
account the anisotropy introduced by the magnetic field 
(40| and have assumed a spherical configuration. Com- 
pact stars with poloidal magnetic fields were studied in 

One can see that the results for the masses and radii 
of the maximum mass stable configuration obtained with 
the inclusion of the magnetic field equal to 2 x 10^^ G 
and 3 x 10^^ G in the SU(2) version of the NJL model is 
larger than the corresponding SU(2) magnetic free star. 
The radius is still quite small, below 8 km, much smaller 
than the average neutron star value, 12—15 km but the 
mass is larger than the value 1.4 Mq which many neutron 
stars have. As referred before, the gravitational mass for 
the two cases including the magnetic field is larger than 
the baryonic mass due to the inclusion of the B^ /2 term 
in the EoS. 



type 


B, 




Mb max 


R 


eo 


B 




Gauss 


(Me) 


(Mq) 


(km) 


(fm-*) 


(10^3 G) 


NJL SU(3) [11] 





1.47 


1.56 


9.02 


7.52 





NJL SU(2) 





1.29 


1.25 


7.11 


13.5 





(NJL SU(2))b 


2 X 10^^ 


1.69 


1.55 


7.63 


12.0 


0.62 


(NJL SU(2))fl 


3 X 10^^ 


1.80 


1.60 


7.97 


10.8 


0.63 



TABLE II: Quark star properties for the EoSs described in 
the text. 



VII. CONCLUSIONS 

The role of the magnetic field effects on quark matter 
has an importance of its own and it has already been 
exploited in the literature. One of the conclusions is 
that the consequences of Landau quantization are not 
negligible for large values of the magnetic field. Quark 
matter is generally described by the MIT bag model. In 
this work we tackle the inclusion of the magnetic field in 
quark matter described by the SU(2) version of the NJL 
model, known to have more realist features. 

Using the standard large- A^c technique we have evalu- 
ated the effective potential for the SU(2) version of the 
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FIG. 6: a) EoS for quark matter in /3-equilibrium and b) 
resulting stellar properties for the SU(2) version of the NJL 
model both with and without magnetic field and the SU(3) 
version without magnetic field. In the figures, Bo = 10^^ G. 

NJL model in 3+1 dimensions which has allowed us to ob- 
tain its thermodynamical potential at finite density and 
in the presence of an external magnetic field. Then, we 
have reviewed some of the most important well known 
results related to three different situations: a) ^ 
and B — (chiral symmetry breaking/restoration), b) 
fi = and B ^ (magnetic catalysis) and c) fi and 
B ^ (magnetic induced quark matter stability). The 
quantity E/A^ for magnetized quark matter described by 
NJL SU(2) model, has a minimum which is lower than 
the one determined for magnetic free quark matter. We 
have also obtained that a magnetic field of the order of 
0.2 X 10^^ G barely affects the effective mass as compared 
with the results for matter not subjected to the magnetic 
field and for i? = 5 x 10^^ G matter is totally polarized 
for chemical potentials below 490 MeV. For small values 
of the magnetic fields the number of filled LL is large and 
the quantization effects are washed out, while for large 
magnetic fields the chiral symmetry restoration occurs 
for smaller values of the chemical potentials. 

In order to introduce /3-equilibrium, we have extended 
the thermodynamical potential so as to consider leptonic 
contributions. Our numerical results show that, for the 
SU(2) case, only very high magnetic fields {B > 10^^ G) 
affect the EoS in a noticeable way. 

While for quark matter in /3-equilibrium, the densities 
at which the pressure is zero increase with the inclusion 
of the magnetic field, the opposite happens in matter not 



subjected to /3-equilibrium. 

Regarding the TOV equation results, we can see that 
the SU(2) version with magnetic field provides masses 
and radii for the maximum mass stable configuration 
larger than the analogous magnetic free configuration. 
The radii are quite lower than the corresponding one 
obtained within the magnetic free SU(3) NJL. This is 
mainly due to the softness of the SU(2) EoS at the lower 
densities. The mass-radius results were obtained with a 
density dependent magnetic field, which increases from 
the surface to the interior of the star. The masses of the 
maximum mass stable configuration for fields of the order 
or larger than 10^^ G are larger than 1.4 Mq, and much 
larger than the analogous magnetic free configuration, 
due to the contribution of the magnetic field. 

The inclusion of the magnetic field in the SU(3) version 
of the NJL is not trivial and is currently under way. We 
expect to obtain larger results for the radius. Knowing 
that the effects of the anomalous magnetic moments is 
very relevant we also intend to take them into account 
in the next calculations. An alternative model, usually 
called PNJLjill , can also be investigated in the presence 
of magnetic fields. 
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APPENDIX A: SUMMING MATSUBARA 
FREQUENCIES AND RELATED FORMULAS 

In this appendix we give the results for the main in- 
tegrals and Matsubara sums appearing along the text. 
The Matsubara sums which are relevant for the different 
inte gral s considered in our work can be derived as (see 
e.g. 0): 



= Ep + Tln{l + e-I^^+fl/^} + Tln{l + e-I^"-^!/"^} , 

(Al) 

which, as T ^ 0, becomes 

lim T ^ ln[(^„ - ZAi)' + 
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+ [/i - Ep]e{^l - Ep) = maxiEp, (A2) 



One can now get rid of the above divergence by adding 
and subtracting a vacuum term (see Eq. (FF)) 



APPENDIX B: EVALUATION OF DIVERGENT 
INTEGRALS 

Consider the divergent term 

^^(2-<5fco)(k/|i?) / ji^MB), (Bl) 

f=uk=0 °° ^ ^ 

appearing in Eq. (fT4)) . By adding and subtracting a LLL 
term to it one gets 



27r 



The subtracted term can be conveniently treated by per- 
forming a change of variables [p']'^ / {2qB) and 
M2 x = ]\/P/{2qB). Then replacing Nf by T,%u 
and performing the integration in = 3 — e dimensions 
leads to 



d oo 



f=u k=0 



dpz 

(27r) 



(B2) 

Now, the integrals can be performed in the following way. 
We change dimensions from 1 ^ d = 1 — e use the stan- 
dard dimensional regularization formula [43| 



27r' 



f=u 



2 2 2 2 

- + Y(l-7i^)-ylnx+- 



(B8) 

Finally, adding all contributions one may write Eq. (|Bip 
as 



oo ,d 



(Tq 



[q' + M'] 



T[A - d/2] 



(47r)''/2r[A](Af2)^-'i/2' 

(B3) 

obtaining, after defining x — AP/{2qB) 



2n'^ r[-l + 6/2] 



/=u fc=0 
1 



(4t 



L/2] 



[A: + a;]-i+^/2 2a;-i+^/2 



(B4) 



Then, using the definition of the Riemann-Hurwitz zeta 
function one can get rid of the summation over Landau 
levels obtaining 



N 

f=u 



a;) In X - 



(B9) 

where the added quantity, P'""-'^^ can be evaluated using a 
sharp cut off reproducing the results quoted in the text, 
which are in agreement with Ref. [2^. Also, note that 
our strategy avoids possible regularization complications 
introduced if one first performs a proper time type of 
calculation and then introduces a three momentum cut 
off. 

The following relation has also been used in the above 
derivation 13111 



0E(l'^/l^)'r[-i+./2] 



1 



C(-l + e/2,x)--^ 



(B5) 

After expanding around 6 = and canceling few terms 
one has the (still divergent) contribution 



^- + \{l-lE)-^\nx-C{-l,x) 

(B6) 



C(-l,a^) 



B2{x) 



x + x^ 



(BIO) 



where Bn{x) represents the Bernoulli polynomial. Fur- 
ther, to obtain the gap equation for the B ^ case 
the following relations 31] are useful, dC,[0,x]/dz — 
lnr(a:)-(l/2)ln(27r), C[o\x] = 1/2-a;, and dC[z,x]/dx = 
— + 1,2;]. 
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